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Abstract 



The relativistic six-quark equations are found in the framework of the dispersion relation tech- 
(Nl ■ nique. The approximate solutions of these equations using the method based on the extraction 
?h ! of leading singularities of the amplitudes are obtained. The relativistic six-quark amplitudes of 
hexaquarks including the quarks of three flavors (u, d, s) are calculated. The poles of these 
amplitudes determine the masses of six-quark systems. 
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In 1977, Jaffe [1] studied the color-magnetic interaction of the one-gluon-exchange potential in 
the multiquark system and found that the most attractive channel is the flavor singlet with quark 
content u 2 d 2 s 2 . The same symmetry analysis of the chiral boson exchange potential leads to the 
CN similar result [2]. 

However, the deuteron channel is not a channel with strong attraction in any baryon interaction 
^ model. If the deuteron had not been found experimentally, it seems highly unlikely that any model 
Q \ would have been able to predict it to be a stable dibaryon. 

It is shown [3] that there are three types of baryon-baryon bound states. The states of the 
>• ■ first type are called deuteron-like states. If chiral fields can provide enough attraction between 
^ , interacting baryons, the systems would be weakly bound. The states of the second type such as 
AA, S*A are named as AA-like states. Due to highly symmetric character in orbital space, these 
systems could be relatively deeply bound, but the strong decay modes of composed baryons cause 
the width of the states much broader. 

The states of the third type are entitled as f2f2-like states. Due to the same symmetric character 
shown in the systems of the second type and the only weak decay mode of composed baryons, 
for instance flfl, these states are deeply bound states with narrow widths. The states of latter 
two types are most interesting new dibaryon states and should be carefully investigated both 
theoretically and experimentally [4-8]. 

There were number of theoretical predictions by using various models [3], the quark cluster 
model [10], the quark-delocation model [11, 12], the chiral SU(3) quark model [13], the flavor 
SU (3) skyrmion model [14] . Lomon predicted a deuteron-like dibaryon resonance using R- matrix 
theory [15]. By employing the chiral SU(3) quark model Zhang and Yu studied flfl and T>Q states 
[16, 17]. 

In the series of papers [18 - 22], a practical treatment of relativistic three-hadron systems 
has been developed. The physics of the three-hadron system is usefully described in terms of 
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the pairwise interactions among the three particles. The theory is based on the two principles 
of unitarity and analyticity, as applied to the two-body subenergy channels. The linear integral 
equations in a single variable are obtained for the isobar amplitudes. Instead of the quadrature 
methods of obtaining solution the set of suitable functions are identified and used us a basis for 
the expansion of the desired solutions. By this means the coupled integral equations are solved 
in terms of simple algebra. In the recent papers [23 - 25], the relativistic three-quark equations 
for the excited baryons are found in the framework of the dispersion relations technique. We have 
used the orbital-spin-flavor functions for the contribution of integral equations. We searched for 
the approximate solution of integral three-quark equations by taking into account two-particle and 
triangle singularities, all the weaker ones being neglected. If we considered such an approximation, 
which corresponds to taking into account two-body and triangle singularities, and defined all the 
smooth functions in the middle point the physical region of the Dalitz plot, then the problem 
was reduced to solving a system of simple algebraic equations. We calculated the mass spectra of 
excited baryons using the input four-fermion interaction with the quantum numbers of gluon [26]. 

In the present paper the relativistic six-quark equations are found in the framework of coupled- 
channel formalism. The dynamical mixing between the subamplitudes of hexaquark are consid- 
ered. The six-quark amplitudes of hexaquarks are calculated. In Sec. II the relativistic six-quark 
equations are constructed in the form of the dispersion relation over the two-body subenergy. The 
approximate solutions of these equations using the method based on the extraction of leading 
singularities of the amplitude are obtained. Sec. Ill is devoted to the calculation results for the 
hexaquark mass spectra (Tables I and II). In the conclusion, the status of the considered model 
is discussed. 

II. Six-quark amplitudes of the hexaquarks. 

We derive the relativistic six-quark equations in the framework of the dispersion relation tech- 
nique. We use only planar diagrams; the other diagrams due to the rules of 1/N C expansion [27 
- 29] are neglected. The current generates a six-quark system. The correct equations for the 
amplitude are obtained by taking into account all possible subamplitudes. It corresponds to the 
division of complete system into subsystems with a smaller number of particles. Then one should 
represent a six-particle amplitude as a sum of 15 subamplitudes: 

A = J2 A ij ■ (1) 

i<3 
i=l 

This defines the division of the diagrams into groups according to the certain pair interaction 
of particles. The total amplitude can be represented graphically as a sum of diagrams. We need 
to consider only one group of diagrams and the amplitude corresponding to them, for example 
A 12 . We shall consider the derivation of the relativistic generalization of the Faddeev-Yakubovsky 
approach. 

In our case the S-wave hexaquarks are considered. We take into account the pairwise interaction 
of all six quarks in the hexaquark. 

For instance, we consider the l uu -diquarks with spin-parity J p = 1 + for the hexaquark con- 
tent (uuuuuu) (Fig. 1). The set of diagrams associated with the amplitude A± 2 can further 
be broken down into three groups corresponding to subamplitudes: A\ m (s, S12345, S1234, S123, S12), 

a 1 UU 1 UU / \ A 1 uu 1 UU-] UU / \ TT 1 1 1 

A L 2 (s, Si2345, S1234, «i2, S34), ^3 {s, §12345, «i2, s 3 4, s 56 ). Here s ik is the two-particle suben- 

ergy squared, corresponds to the energy squared of particles i, j, k, s^i is the energy squared 
of particles i, j, k, /, corresponds to the energy squared of particles i, j, k, I, m and s is the 

system total energy squared. 

The system of graphical equations is determined by the subamplitudes using the self-consistent 
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method. The coefficients are determined by the permutation of quarks [30, 31]. In order to 
represent the subamplitudes Ai, I = 1 — 3 in the form of a dispersion relation, it is necessary 
to define the amplitudes of quark-quark interaction. The pair quarks amplitudes gg — > qq are 
calculated in the framework of the dispersion N/ D method with the input four-fermion interaction 
[32 - 34] with the quantum numbers of the gluon [26, 35]. The regularization of the dispersion 
integral for the D-function is carried out with the cutoff parameter A. 
The four-quark interaction is considered as an input: 



g v (qXIf^qf + 2 gfi (qXIf^q) (sX-f^s) + g { y S) (sX-y^sf . (2) 

Here If is the unity matrix in the flavor space (u, d), X are the color Gell-Mann matrices. Dimen- 
sional constants of the four-fermion interaction gy, gy\ and gy^ are parameters of the model. 

At gy = gft = gy^ the flavor SU (3) f symmetry occurs. The strange quark violates the flavor 
SU (3) f symmetry. In order to avoid additional violation parameters we introduce the scale of the 
dimensional parameters [35]: 

A= mk) 

(rrii + m k y 

Here rrii and m k are the quark masses in the intermediate state of the quark loop. Dimensionless 
parameters g and A are supposed to be constants which are independent of the quark interaction 
type. The applicability of Eq. (2) is verified by the success of De Rujula-Georgi-Glashow quark 
model [26], where only the short-range part of Breit potential connected with the gluon exchange 
is responsible for the mass splitting in hadron multiplets. We use the results of our relativistic 
quark model [35] and write down the pair quark amplitudes in the form: 

«n(«») = ; ° 2 f f , , (5) 

1 — r>n{Sik) 



(m^+m fc ) 2 A 



Bj Slk )= [ ^ih P^y^ik) (6) 

, J , a S lk -S ik 

Here G n (sik) are the diquark vertex functions (Table III). The vertex functions are determined 
by the contribution of the crossing channels. The vertex functions satisfy the Fierz relations. 
All of these vertex functions are generated from g v , g^ and gy S \ B n (sik) and p n (s ik ) are the 
Chew-Mandelstam functions with cutoff A [36] and the phase spaces, respectively: 



Pn(s ik ,J PC ) = (a(n,J PC ) . g - + fj(n, J PC ) + 8(n, J PC ) K ^ 
\ (mi + m k y s ik 



v {sik - (m + m k ) 2 )(s ik - (m 8 - m k ) 2 ) 

&ik 

The coefficients a(n, J PC ), (3(n, J PC ) and 5(n, J PC ) are given in Table III. 
Here n — 1 coresponds to gg-pairs with J p = + , n = 2 corresponds to gg-pairs with J p = 1 + . 
In the case in question the interacting quarks do not produce a bound state, therefore the 
integration in Eqs. (8) - (10) is carried out from the threshold (rrii + ^k) 2 to the cutoff A(ik). 
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We consider the hexaquark state with the strangeness S = 0, the isospin 1 = 3 (spin-parity 
J p = 0+2+). 
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The coupled integral equations correspond to Fig. 1 can be described as: 



A{ UU {S, S12345, S1234, S123, S12) = r/^D ^/ 12 \i + 8 ^l( S 12, 1 UU )A\ UU (S, S12345, S1234, Sl23, S' l3 ) 

[I — rf 1™ [s 12 )\ 



A 2 (s, S12345, S1234, S12, S34) 



. ^ 14 LI LL U, ^ Li LL , 

^3 (S, 512345) S12, S34, S56J = 



, •S123455 s 1234, S13, S24) ) (8) 
\2B1uu yS\2)B\uu ( s 34) 

[1 - B 1 uu(s 12 )][l - Si^(s 34 )] 
+ 4J 3 (s 12 , S34, 1" M , (s, S12345, S1234, si 23! 4s) 

+ 8 J4(Si2, S34, 1™, l" M )^i ( s ) s 12345) s 12355 s 125, s 'l5) 

i 5 12356> 5 1256; S 15! S 2d) 
+ 8J 6 ( Sl2 ,S34,l Un ,l UU )^ 1 (S 1 Sl2456? s U56i s l5i s 4q} 
+ 16J 7 ( Sl2 ,S34,l n ",l n ")^2 i s 12345; s 2345; s 23' ^45) 

+ 8J 8 ( Sl2 ,S34,l UM ,l UM )^3 ; 5 12345; s 15? s 23' S 46) ' (^) 

X 3 Biuu {Si2)Biuu{S34)Biuu (ssq) 



where 



[1 - B 1 uu(si 2 )][l - B 1 uu(s 34 )][l - Bi««(s 56 )] 

+ 12J 9 (S 12 , S34, S 56 , 1™, I"", l W )Ar>, S12345, S1234, S'i2 3 , 4 3 ) 

+ 24 Ji (si2 ) s 34, S56, J- , 1 j s 12345; s 2345; s 23> S 45/ > 

(10) 



(m 1 +m 2 ) 2 A i 

.x Gi(5i 2 ) / ds' 12 Gi(s' 12 )p i (s' 12 ) + r dzijl) / in 

[1 - Biisu)] J tc s' 12 - S12 J 2 

(mi+m2j — 1 

(m 1 +m 2 ) 2 A i ^ ^ 

Gi(si 2 ) / ds' 12 G i( s 12) Pi( s 12) 1 /* dz i( 2 ) f dz 2 (2) 

J2{Sl2,t) = 77 7J7 77 / ; — 

[1 - Bi(s 12 )\ J 71 s' 12 - S12 2tt 

(mi+m2j — 1 — 1 

^3(2)+ 

x / dz 3 (2) . (12) 

J 2) _ y/1 - z\{2) - z\(2) - 4(2) + 2z 1 (2)z 2 (2)z 3 (2) 
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(m 1 +m2) Aj 

w . .x ^(312)^(^34) / ds' 12 Gi{s' l2 )pi(s' l2 ) 

J3(Sl2 ' S34 ' 2 ' Jj " [1 - B t (s 12 )}[l - Bj (s u )} J n s' l2 -s 12 " 

[m\+m2) 

(m 3 +m 4 ) 2 Aj 

x / ^ 4 G , J (s / 34 )p J (s / 34 ) + r 1 dz 1 (3) J dz 2 (3) 
J it S34 — S34 J 2 J 2 

(m 3 +m4) 2 —1 -1 

f / ■ -\ -^i( s 34) i / • \ 

J5{Sl2,S U ,t,3) = _ - -^-^ J 2 (Sl2,») , 

Je(si2,S34,i, j) = il(si2,«) • A(S34, j) , 

f/ „ • ■ \ _ gi(si2)Gj(s34) / ds' 12 Gi(s' 12 )pi(s' 12 ) 

H l2, 34 " j) " [1 -%)][!- ^.( S34 )] , 7 S ' 12 - Sl2 

(mi+m2)' 1 

(m 3 +m 4 ) 2 Aj 

/ ^34 G i(4)Pi(4) 1 /'^i(7) /-^ 2 (7) /" 6^(7) 



J it sL - s 34 2vr 7 2 7 2 7 2 

(m 3 +m 4 ) 2 -1 -1 -1 



2 4 (7) + 

^1 - zf(7) - 4(7) - zj(7) + 2z 1 (7)z 3 (7)z 4 (7) 



x / ^ (7) 7T^ 



24 (7) 



(m 1 +m 2 ) A, 

G i(s 12) Gj(s 34) f ds' 12 Gi(s' l2 )pi(s' 12 ) 



J, ■ -\ _ <J"i( 1 Sl2j<- T A' S 34j /" US 12 Crj(S- 

8lSl2 ' S34 ''' jj " [l-^( Sl2 )][l-^( S3 4)] , 7 ^ 1? ^ 2 - Sl2 



(m 3 +m 4 ) 2 A J - 



X 



/■ ds. 



(is 3 4 ( ^j(' S 34)Pi(' S 34) 



(m 3 +m 4 ) 2 34 

+1 +1 +1 24(8)+ +1 z 6 (8)+ 

x JL /^iM /^M /^8) I ^ (g) f dz^ f ^ 



(2vr) 2 _{ 



2J 2 J 2 J J 2 

-1 -1 24(8)- -1 2 6 (8) 



X 



X 



^/l - ^(8) - 5(8) - ^(8) + 2^(8)^3(8)^4(8) 

1 

yjl - z 2 2 (8) - 3(8) - 5(8) + 2^(8)Z5(8)^ 6 (8) 



(Sl2,S34,S56,«,J, «) = Tj 75-7 77^3 («12, S34, I, j) , 

[1 - B k {S^)\ 



(m 1 +m 2 ) 2 A, 

G f »(si2)G ! j(s34)G f fc (s5 6 ) r ds' 12 G , i (s' 12 )p i (si 2 ) 



i , . . , \ ^i{Sl2)^j{S34)^k{S56) f <tS' v 

Jw{si2,su,sm,i,J,K) — t, 57 — 77 57 — 77. 57 — 7 / , 

[1 - £i(Si 2 )][l - 5 J -(S 34 )][1 - 5jfc(s 5 6j] y 7T s' 12 - S12 

(mi+ni2) 

(m 3 +m 4 ) 2 Aj (m 5 +m 6 ) 2 A fc 

x / ^g3igj(£34)M£34) / ^56 <^fc(46)Pfc( g 56) 

J TV S' u - S34 J TV S' 56 - S 56 

x 1 J dz.jlO) J dz 2 (l0) J dz 3 (10) J dz A (10) Z5 J )+ dz5 (W) 

^-l 2 -1 2 "I 2 "I 2 25(1-)- 

x 1 (20) 

^/l - ^ 2 (10) - z 4 2 (10) - zi(10) + 2^(10)^(10)^(10) 

We should discuss the coefficients multiplying of the diagrams in the equations of Fig. 1. For 
example, we consider the first subamplitude Ai(s, S12345, S1234, S123, ^12)- In the Eq. (8) (Fig. 1) 
the first coefficient is equal to 8, that the number 8 = 2 (permutation particles 1 and 2) x4 (we 
can use third, 4-th, 5-th, 6-th particles); the second coefficient equal to 12 that the number 12 = 4 
(used third, 4-th, 5-th, 6-th particles) x3 (in this case we can consider 4-th, 5-th, 6-th particles). 
The similar approach allows us to take into account the coefficients in the Eqs. (9) and (10). 

Let us extract two- and three-particle singularities in the amplitudes A\ uu (s, S12345, S1234, S123, S12), 

A 2 (S, S12345, S1234, Sl2, S34J, /I3 [S, S12345) S12, S34, Ssq): 

A l uu / \ a l ( S 5 s 12345 5 S1234, Sl23) Si2)Biuu(,Si2) , , 

A 1 (S,S12345,S1234,S123,S12) = 71 5 7 yi , (21) 

[1 — Diuu{Si2)\ 

Yunnan < . Cl\ 1 (s, Sl2345 ) ^1234 ; S 12 ; ^34) B\uu (s 12 ) B\uu ( S34) 

a 2 ys, S12345, S1234, S12, S34J - — — — — — , yzz) 

[1 - ±>i««(Si2j [1 - -Dl"H s 34jJ 



A3 [S, Sl2345i s 12) S34, S56J = 



Q; 3 (s, Si2345 ? s 12, S34, S^)B\uu (Si2)^l«« (s34)-Bl«« (S56) 



[1 - Si««(si 2 )][l - -Bi« u (s 34 )][l - -Bi^(s 56 )] 



(23) 



We do not extract four-particles singularities, because they are weaker than two- and three- 
particle singularities. 

We used the classification of singularities, which was proposed in paper [37]. The construction of 
the approximate solution of Eqs. (21) - (23) is based on the extraction of the leading singularities 
of the amplitudes. The main singularities in s ifc = (m^ + m k ) 2 are from pair rescattering of the 
particles % and k. First of all there are threshold square-root singularities. Also possible are pole 
singularities which correspond to the bound states. The diagrams of Fig. 1 apart two-particle 
singularities have triangular singularities and the singularities defining the interactions of four, 
five and six particles. Such classification allows us to search the corresponding solution of Eqs. 
(8) - (10) by taking into account some definite number of leading singularities and neglecting all 
the weaker ones. We consider the approximation which defines two-particle, triangle and four-, 
five- and six-particle singularities. The contribution of two-particle and triangle singularities are 
more important, but we must take into account also the other singularities. 

The functions ai, I — 1 — 3 are the smooth functions of s^, s^, s^i Sijkim as compared with 
the singular part of the amplitudes, hence they can be expanded in a series should be employed 
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further. Using this classification, one defines the reduced amplitudes a>i, a 2 , a 3 as well as the 
-B-functions in the middle point of physical region of Dalitz-plot at the point s : 



6 

s + 4 J2 mf 



hk 



i,k=l 
i<k 



3 3 



S123 = S m2 ik ~ 2 m i > ( 25 ) 



,fc=i i=l 
;<fc 



4 4 

S1234 = SO m lk ~ 2 H m i • ( 26 ) 



i,fc=l 1=1 



Such choice of point s allows us to replace integral equations (8) - (10) (Fig. 1) by the algebraic 
equations (27) - (29), respectively: 

a\ m = X + 8h(r u r u )al uu + 12I 2 (l uu l uu r u )al uuluu , (27) 
af 1 ™ = A + 4J 3 (l uu l u?t l uu )ar + 8/ 4 (l»l m l u >r + u^r^l^a 1 ^ 1 ™ 



+ 8/ 6 (l MU l M "l" M l MU )a^ 1111 ' + 16/ 7 (l" M l uu l MU l" M )a^ luu + 8/ 8 (l« M r«l« M l« M r«)ai ulllllt,lu " , 

(28) 

a l uu ^ uu = a + i2/ 9 (r M i Mu r M r u )ar + 24i w (i uu i uu r u r u r u )ar iuu , (29) 

where Aj are the current constants. We used the functions I±, I 2 , h, h, h, h, I7, h, h, ho- 



(m 1 +m 2 ) 2 A i 



T( -, Bjjsl 3 ) f ds' l2 G^)p,{s' 12 ) f d Zl (l) 1 

1[3) " B^sf) J vr s> 12 -stf J 2 l-^Vis)' 1 J 

(mi+m,2) — 1 



+1 , +1 



_ ^pffwff) / ^ gKfgWsy 1 ?^i(2) 7^2(2) 

MW ~ " B^sf) J vr s' l2 -sf 2vri 2 7 

(mi+rri2) 2 —1 —1 



2 



^3(2)+ 

x / dz 3 {2) 

J 2) - \A - *?(2) - z 2 2 (2) - 5(2) + 2^(2)^(2)^(2) 

x n ETTTTT i ETTTTT ' ( 31 ) 



/ 3 (u'fc) 



5,(^3) f ds' 12 Gj{sf) Pi {s' 12 ) 



(mi+m2) 



8 



I 8 (ijklm) 



(m 3 +m 4 ) 2 Aj 



X 



ds', A G](sf) Pj (s^) + rd Zl (3) ) dz 2 (3) 



(m3+m,4) 2 -1 —1 



J 7T Soa Sr^ J 



s'^-sf 7 2 7 2 1-B fc (4)' 

h{ijk) = h(ik), 
h(ijkl) = I 2 {ikl) , 
hiijkl) = h{ik) ■ h(jl) , 



I r (ijkl) = 



x 



X 



X 



X 



X 



X 



(m 1 +m 2 ) 

B^B^sf) ' } ds' 12 G^) Pi (s' 12 ) 



(m 3 +m 4 ) 2 A ■ 



4^s 34 2tt7 2 7 2 7 2 



ds' u G)(sf) Pj (s' u ) 1 + rd Zl (7) fdz 2 (7) fdz 3 (7) 



7 7T Sq 4 — S34 27T 7 



(m 3 +m 4 ) 2 -1 -1 -1 

z 4 (7) + 



J dz,(7)j= 



J } - \A - *?( 7 ) " W) - 4(7) + 2^(7)^(7)^(7) 

1 1 

i - ^(sy i - s^y ' 

(m 1 +m 2 ) 2 A i 

^(^ 5 )^(^ 3 )S m (^ 6 ) """ f ds' l2 G*(sf) Pl (s> l2 ) 



Bi{sl 2 )Bj(sf) 7 7T S ' 12 -^ 2 

(mi+m 2 ) 2 



(m3+m4) 2 Aj 



^ g(fjj%(fk) 
(m 3 +m 4 ) 2 

+ 1 , /ox +1 , /^x +1 , /on 2 4(8)+ +1 2 6 (8)+ 



7 7T 



-1 -1 -1 24(8)- -1 z 6 (8)- 



1 - zf(8) - z 2 3 (8) - zl(8) + 2zi(8)z 3 (8)z 4 (8) 
1 

y/l - zj(8) ~ 4(8) - 4(8) + 2z 2 (8)z 5 (8)^ 6 (8) 
1 1 1 



l-^yi-^^l-S^)' 
I 9 (ijkl) = h(ijl) , 



(m 1 +m 2 ) 2 A i 



BtfflB^tf) f ds' 12 G^) Pi {s' l2 ) 



I „12 



X 



B^B^B^sf) J 7T s 12 -o 

(m\+m2) 

(m 3 +m 4 ) 2 Aj (m 5 +mg) 2 A fc 

I ds' u G%sf) P] {sy) I ds> 5e Gl(s™)p k (J 56 ) 

J 71 S34 Sq J 71 Sgg Sq 

(m3+m4) 2 (m 5 +m6) 2 

x 1 J dgi(io) J ^ 2 (io) J ^3(io) J ^4(10) Z5 J )+ dZ5{10) 
n -1 -1 -1 -1 



«6(1-)- 



X 



x n p~77TT i D w (39) 



^1 _ 2 2(io) - 5(10) - z|(10) + 2^(10)^(10)^5(10) 
1 1 



where i, j, k, I, m correspond to the diquarks with the spin-parity J p = + , 1 + . 

In the equation (30) 21 (1) is the cosine of the angle between the relative momentum of particles 

1 and 2 in the intermediate state and the momentum of the particle 3 in the final state taken in 
the cm. of particles 1 and 2. We can go from the integration of the cosine of the angle dzi(l) to 
the integration over the subenergy ds' l3 . 

In Eq. (31) 21 (2) is the cosine of the angle between the relative momentum of particles 1 and 

2 in the intermediate state and the momentum of the particle 3 in the final state taken in the 
cm. of particles 1 and 2, z 2 (2) is the cosine of the angle between the momenta of particles 3 and 
4 in the final state of cm. of particles 1 and 2, -23(2) is cosine of the angle between the relative 
momentum of particles 1 and 2 in the intermediate state and the momentum of the particle 4 in 
the final state of cm. of particles 1 and 2. Then we pass from dzi(2)dz2(2)dzs(2) to ds' 13 ds' u ds' 24 . 

In Eq. (32) 21 (3) is the cosine of the angle between the relative momentum of particles 1, 2 
in the intermediate state and the relative momentum of particles 3, 4 in the intermediate state 
in cm. of particles 3 and 4; 2 2 (3) is the cosine of the angle between momentum of particle 3 in 
the intermediate state and relative momentum of particles 1, 2 in the intermediate state in cm. 
1 and 2. We pass from dz 1 (3)dz 2 (S) to ds' 123 ds' 23 . The similar method are used for the functions 
(33) - (35), (38). 

In Eq. (36) 21 (7) is cosine of the angle between relative momentum of the particles 1, 2 in the 
intermediate state and the relative momentum of particles 3, 4 in the intermediate state in cm. 
of particles 3 and 4; ^(7) is the cosine of the angle between the momentum of particle 3 in the 
intermediate state and relative momentum of particles 1, 2 in the intermediate state in cm. of 
particles 1 and 2; z 3 (7) is cosine of the angle between momentum of particle 5 in the final state and 
relative momentum of particles 1, 2 in the intermediate state in cm. of particles 3 and 4; -24(7) is 
cosine of the angle between momentum of particle 5 in the final state and the relative momentum 
of particles 3, 4 in the intermediate state in cm. of particles 3 and 4. Then we translated the 
dzi(7)dz 2 (7)dz 3 (7)dz4 : (7) to ds' 123 ds' 23 ds' 125 ds' 45 . 

In Eq. (37) 21(8) is the cosine of the angle between momentum of particle 5 in the final 
state and the relative momentum of particles 1, 2 in the intermediate state in cm. of particles 
1 and 2; 2:2(8) is the cosine of the angle between the relative momentum of particles 1, 2 in 
the intermediate state and the relative momentum of particles 3, 4 in the intermediate state in 
cm. of particles 3 and 4; 2:3(8) is the cosine of the angle between momentum of particle 3 in 
the intermediate state and the momentum of particle 5 in the final state in cm. of particles 1 
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and 2; 2:4(8) is the cosine of the angle between the momentum of particle 3 in the intermediate 
state and the relative momentum of particles 1, 2 in the intermediate state in cm. of particles 

1 and 2; 25(8) is the cosine of angle between momentum of particle 6 in the final state and the 
relative momentum of particles 1, 2 in the intermediate state in cm. of particles 3 and 4; 26(8) 
is the cosine of the angle between momentum of particle 6 in the final state and the relative 
momentum of particles 3, 4 in the intermediate state in cm. of particles 3 and 4. We pass from 
dzi(8)dz2 (8)^23(8)^24(8)^25(8)^26(8) to ds\ b ds' 12 ?,ds'2 b ds'2zds' 12& ds' m . 

In Eq. (39) 21 (10) is the cosine of angle between relative momentum of particles 1, 2 in the 
intermediate state and the relative momentum of particles 3, 4 in the intermediate state in cm. 
of particles 3 and 4; 22(10) is the cosine of angle between the relative momentum of particles 1, 

2 in the intermediate state and momentum of particle 3 in the final state in cm. of particles 1 
and 2; z 3 (10) is the cosine of the angle between the relative momentum of the particles 3, 4 in the 
intermediate state and the relative momentum of particles 5, 6 in the intermediate state in cm. 
of particles 5 and 6; 2:4 (10) is the cosine of angle between relative momentum of particles 1, 2 in 
the intermediate state and the momentum of particle 5 in the final state in cm. of particles 3 
and 4; 25(10) is the cosine of the angle between the relative momentum of the particles 3, 4 in the 
intermediate state and the momentum of particle 5 in the final state in cm. of particles 3 and 4. 
We pass from dz 1 (10)dz 2 (10)dz 3 (10)dz 4: (10)dz 5 (l0) to ds'^ds^ds'^ds'^ds'^. 

The other choices of point Sq do not change essentially the contributions of ai, I = 1—3, 
therefore we omit the indices s k . Since the vertex functions depend only slightly on energy, it is 
possible to treat them as constants in our approximation. 

The solutions of the system of equations are considered as: 

a i (s) = F i (s,\ i )/D(s), (40) 

where zeros of D(s) determinants define the masses of bound states of dibaryons. 

As example, we consider the equations for the quark content uuuuuu with the strangeness 
5 = 0, the isospin 1 = 3 and the spin-parity J p = + ,2 + (Fig. 1). The similar equations have 
been calculated for the strangeness 5 = 0, —1, —2, —3, —4, —5, —6, the isospin I — 0, \, 1, §, 2, |, 

3 and the spin-parity J p = + , 1 + , 2 + . We take into account the u, d, s quarks. 

In Appendix I the NN S u=ooi, AA S /j =0 oi, tt£lsiJ=-mo, ^\su=-2oo, N ^su=-3±2 dibaryons are 
given. 

III. Calculation results. 

The poles of the reduced amplitudes (i — 1 — 3) correspond to the bound state and determine 
the mass of the hexaquark with the quark content (uuuuuu), with the isospin 1 = 3 and the spin- 
parity J p = + ,2 + . The quark masses of model m U) d = 410 MeV and m s = 557 MeV coincide 
with the ordinary baryon ones in our model [38]. 

The model in question has only three parameters: the cutoff parameter A = 11 (similar to 
the model [41]) and the gluon coupling constants g and g 1 . These parameters are determined 
by the AA and the di-Q masses. We have considered the two type of calculations. In the first 
case we use the gluon coupling constants g\ = 0.292 (diquark 1 + ) and go = 0.653 (diquark + ). 
Which are fitted by the AA state with the M = 2173 MeV and the di-ft with the M = 3232 MeV , 
respectively. In the second case the gluon coupling constants g\ = 0.325 and go = 0.647 are 
determined by the masses of AA state with the M = 2171 MeV and the di-ft state M = 3093 MeV. 
The experimental data of these masses are absent, therefore we use the paper [12]. In our model 
the correlation of gluon coupling constants g and g\ is similar to the S-wave baryon ones [38]. 

The estimation of theoretical error on the S-wave hexaquarks masses is 1 MeV . This results 
was obtained by the choice of model parameters. We predict the deuteron state as the mix of S- 
and D-wave contributions. (NN SI j =om with the mass M = 1865 MeV and AA s/ j =001 with the 
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mass M = 1834 MeV). The experimental value of deuteron mass is M = 1876 MeV. In the cases 
of the NN, A A, NA systems the Pauli principle requires that (-1) L+/+J = (-1), where L the 
orbital moment, I isospin, J spin of state are respectively. The wave function of dibaryon must be 
antisymmetric for the permutation of all quarks. If we consider the generalized Pauli rule for the 
wave function of dihyperons, we must add the strangeness contribution to the isospin I + | . Then 
we obtain the formula ( — l) i+/ +f +J = ( — 1). This rule allows us to suggest the classification of 
dibaryons with the certain strangeness, isospin and spin-parity (Tables I and II). We predict the 
degeneracy of the some states. The contributions of subamplitudes to the hexaquark amplitude 
are shown in the Appendix I (for example, AAs7j = _ 2 oo)- The nonstrange dibaryon with the isospin 
1=1 and the spin-parity J p = + is absent. 

The states NA and AA with the isospin 1 = 1 and the spin-parity J p = 2 + possess the mass 
M = 2020 MeV. For the NA and A A with the isospin 1 = 2 and the spin-parity J p = 1 + we 
obtained the mass M = 1984 MeV . For the state AA with the isospin 1 = 3 and the spin-parity 
J p = 0+, 2+ (M = 2379 MeV) the degeneracy is predicted. It is shown in the Table I. 

The results for the strange sector of model are given in Table I and II. 

IV. Conclusion. 

The dibaryon physics can be very delicate [39 - 41]. The deuteron channel is not a channel 
with strong attraction in any baryon interaction model. If the deuteron had not been found 
experimentally, it seems highly unlikely that any model would have been able to predict it to be 
a stable dibaryon. 

The H-particle (SIJ = —200) is a six quark state consisting mainly of octet-baryons, similar 
to the deuteron and one can find only a weak attraction in the model [41]. Hence, a qualitative 
analysis is insufficient to judge whether or not the H-particle is strong interaction stable. System- 
atically, the authors find that a strong attraction develops only in decuplet-decuplet channels and 
a mild attraction in octet-decuplet channels [41]. Moreover, in the H-particle case, the channel 
coupling effect may even be more important than the deuteron case. In fact, it is bound without 
taking coupled channels into account. Besides the binding energy of the H, an interesting question 
regarding the H is its compactness, i.e. whether the H is a compact 6-quark object or a loosely 
bound AA state. 

For systems with strangeness S = —3, Pang et al. have calculated the state NQ(SIJ = —3^2), 
which was shown to be midly attractive, with energy below AEir threshold. 

They have carried out a dynamical channel coupling calculation to examine this state futher. 
The NQ, AS*, HE*, SS* channels are all included. The authors find this to be a compact 

six quark state [41]. 

For systems with S = —4, with the quantum numbers SIJ = —410 as an example, the lowest 
mass channel is composed of two octet baryons from the same isodublet. The result shows that 
the system with S = —4, / = 1, J = is unbound, even when the and E*f2 channel couplings 
are taken into account. 

For comparison, Pang et al. have also calculated the SIJ = —401 state. The SS, SS*, AQ 
and S*S* coupling channels are included. The result is very similar to the SIJ = —410, i.e. they 
do not find a bound state in this channel. 

For the systems with S = —5, Pang et al. take the SIJ = — 5|0 state as an example. This state 
is interesting as a Pauli principle favored state. If only two-baryon S-wave channels are taken into 
account, there is only one channel for this state. The calculation shows [41] that the contribution 
of the kinetic energy term, due to quark exchange and derealization effects, contributes strongly 
towards the formation of a bound state. However, the one-gluon-exchange interaction largely 
compensates for this attraction. Pang et al. conclude that this state is not a good candidate for 
a dibaryon resonance search due to its small binding. 
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In the paper [41], Pang et al. would present a systematic study of possible candidates of S-wave 
baryon-baryon bound states. 

The H-particle, TVQ-state and di-fl may be strong interaction stable. Up to now, these three 
interesting candidates of dibaryons are still not found or confirmed by experiments. It seems that 
one should go beyond these candidates and should search the possible candidates in a wider region, 
especially the systems with multi-strangeness, in terms of a more reliable model. 

In our model the deuteron consist of the A A, NN contributions (the strangeness S = 0, the 
isospin 1 = 0, spin-parity J p = 1 + , the quark content is uuuddd). The H-particle (SIJ = —200) 
content includes NE, E*E*, EE, AA. 

For the systems with strangeness S = -3 (Nfl SIJ = — 3|2) the Nfl, AS*, HE*, ES*, S*E* 
channels are included. 

For the di-fl state we consider the strangeness S = — 6 (SIJ = —600). 

The gluon coupling constants in our model is determined by the masses of the H-particle and 
di-fl state (Table I and II). 

We considered 39 dibaryons, calculated the masses these states with the strangeness S = 
0, -1, -2, -3, -4, -5, -6 and the spin-parity J p = + , 1 + , 2+ 

In our paper the dynamics of quark interactions in defined by the Ghew-Mandelstam functions 
(Table III). We include only three parameters: the cutoff A, gluon coupling constants g , g\. The 
relativistic six-body approach gives rise to the dynamical mixing of the six-quark amplitudes and 
the dibaryon amplitudes. We calculated the masses of two groups of dibaryons (Table I and II), 
which similar to the results of other papers [12, 38 - 40]. In our paper the relativistic description 
of six particles amplitudes of S-wave dibaryons are considered. We use only three parameters for 
the calculations of 39 dibaryon masses. The interesting research is the consideration of the qqqqqQ 
states with Q a heavy quark (Q = c,b). 
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Appendix I. The reduced amplitudes of dibaryons NN S u=ooii ^^su=ooii ^5/j=-600? 



NN 



5/7=001- 



= A + 2a i ri 1 (l uu r u ) + 6a? h(l uu ud ) + 6 I 2 (l uu l uu ud ) + 6 a° 2 l 2 (l uu ud ud ) 

af = A + 2 al <w / 1 (l dd l d(i ) + 6 h(l dd ud ) + 6 a^ d °" d / 2 (l dd l dd ud ) + 6 af QUd l 2 (l dd ud ud ) 

< d = \ + 2arh(0 ud r u ) + 2a l "h(Q ud l dd ) + Aaf h(0 ud ud ) + h(0 ud l uu l dd ) 

+ 2 ar°" d / 2 (0 w r"0" d ) + 2 a 1 2 M ° ud I 2 (0 ud ud l dd ) + 4 a° ud °" d / 2 (0 w w w ) 

a i««i«« = A + 0,0^(4/3(1^1^0^) + 2 h(i uu i dd o ud ) + 2 h(i dd r u o ud )) + 2a 2 uu ° ud i 5 (r u i dd r u o ud ) 



+ 2a 1 2 dd ° ud I 5 (l M r u ud l dd ) + af d °" d (4 J 6 (r u l d W d ) + 4/7(1^1^0^0^) 

+ 4 l 7 (l dd l uu ud ud )) 

al uu ° ud = A + 2a{ uu I 3 (l uu ud l uu ) + 2 a{ M h(0 ud l uu l dd ) + a?" "(2 J 3 (l" M Md Mc( ) + 4 J 4 (l M "0" d 0" d ) 

+ 2 J 4 (0 Md l M "0" d )) + 4 a i uu i dd / 7 (i««o nd l" n l dd ) + 4 a^ ,10Ud 7 7 (0 u<i l uu l uu ud ) 

+ a^°" d (2/ 5 (0 ud l uu u<i l dd ) + 4/ 6 (l uu ud u<i l dd )) + d ° ud (2/5(l MU U(i ud ud ) 

+ 4 / 6 (i««o" d O" d O" d ) + 4 / 7 (i««o ud O nc( Md ) + 4 j 7 (o Md l M "0 Mc( Mc! )) 

+ 4af lWotId / 8 (l™0^0 Md f u l drf ) 

a i dd °" d = A + 2a{ uu h(0 ud l dd l uu ) + 2a{ M h(l dd ud l dd ) + aT" (2 l 3 (l dd ud ud ) + 4 h(l dd ud ud ) 

+ 2 h(0 ud l dd ud )) +4al uuldd I 7 (l dd ud l dd r u ) + al uu ° ud (2 l 5 (0 ud l dd ud l uu ) 

+ AI 6 (l dd ud ud l uu )) +4a 1 2 d ° ud I 7 ( y ud l dd l dd ud ) + a° 2 d ° ud (2 l 5 (l dd ud ud ud ) 

+ 4 / 6 (i rfrf o nd O Md O Mc( ) + 4 i 7 [\ dd Q ud Q ud Q ud ^ _|_ 4 /^O^l^O^O*"*)) 

+ 4 a i u "i dd o" d j 8 ( 1 ^ ^ 1 ^ r «) 

af d °" d = A + af" (7 3 (0 u<i u<i l uu ) + 2 J 4 (0 Md nc 4 MU )) + <"(J 3 (0 ,ld ,ld l <w ) + 2 / 4 (0" d 0" d l dd )) 
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+ a[ ud (2 1 3 (0 ud ud ud ) +4h(O ud O ud O ud )) + {2 1 5 {0 ud ud r u l M ) + 2 I e (O ud O ud l uu l dd ) 

+ 2 i 7 (o ud o ud i dd r u ) + 2 i 7 (o ud o ud r u i dd )) + a l uu ° ud (2 i 6 (o ud o ud i uu o ud ) + 2 i 7 (o ud o ud r u o ud ) 

+ 2 I 7 (G ud a ud a ud l m )) + al dd0U \2 I 6 (G ud wi a ud l dd ) + 2 / 7 ( « + 2 I 7 (G ud a ud l M a ud )) 

+ a° 2 ud ° ud (2 / 5 (o« d O ud O U(i O ud ) + 2 / 6 (o« rf O ud O Md O w ) + 4 J 7 (0 M W d M(i )) 

-|- ^l uu l dd O ud ^ j^^gtirfgwrfj^wMQMrfj^rfd^ _|_ 2 ^^QudQud-^uu-^ciciQUc^ _|_ 2 j^gncfgMfifg'ticf^Mu^cfcJ^ 

a^ ldd °" d = A + 2a{ uu I 9 (r u ud l dd r u ) + 2 a 1 ^ I 9 (l dd ud l uu l dd ) + a° ud (4 l 9 (l uu l dd ud ud ) 

+ 2I 9 (l uu ® ud l dd Q ud ) + 2/9(1^0^10^)) + 4a^ 1 ' ldd / 10 (l UM Mc! l' id l UM l^) 

+ 4 a^° ud / 10 (l dd l uu ud ud l uu ) + 4 ai dd011d / 10 (l M "1^0 uc! M(i l^) 

+ af d ° ud (4/io(r u l dc( O u W d ) + 4/ 10 (i^i««o Md O Md O Md ) + 4 / 10 (l M "0 Md 1^0 Md Mc( )) 

AA 5 /j =001 : 

a?" = A + 2 af/^rr) + 6 af^l^O^) + 6 a r° ud I 2 (l uu l uu ud ) + 6 < d o" d / 2 (i- wd O ud ) 

af = A + 2 af J^l^) + 6 <" J^l^O"") + 6 af " I 2 (l dd l dd ud ) + 6 a^ " h(l dd O ud O ud ) 

aT d = A + 2a{ uu h(0 ud r u ) + 2 h(O ud l dd ) + Aaf h(O ud O ud ) + 4 ag" 1 " W^V") 

+ 2 ar° Ud /2(0 ud l uu u<i ) + 2 af 011d / 2 (0"V d l dd ) + 4 ^^^"^(O^O^O^) 

a i««i«« = A + 2a{ uu h(r u l dd r u ) + 2al dd h(l dd r u l dd ) + aT\4:I 3 (r u l dd ud ) + 2 h(l uu l dd O ud ) 

+ 2 h(i dd r u o ud )) + Aal uuldd i 6 (r u i dd r u i dd ) + al uuoud (2i 5 (r u i M r u o ud ) 

_|_ ^ j^-^dd-^ mtg ltd mi ^ _|_ ^l dd o ud ^ ^cfc(-^«ug«d-^cic^ _|_ ^ j^^uM-^ddgud-^dd^ 

_)_ Q:°" doud (4/ 6 (l" M l< M 0" d O" d ) + 4/ 7 (i« M i (M g« c( g« rf ) -f 4 / 7 (i dd i uu o ud O ud )) 

+ 4 Q; 1 "" lddo " d j g (i« u i^i u * t () ,trf l d<i ) 

al uu ° ud = A + 2a 1 luu / 3 (l" M 0" d l M ") + 2 a^/^O^l""!^) + a?" "(2 I 3 (l uu ud ud ) + 4 J 4 (l M "0" d 0" d ) 
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a 



+ 2 i 4 (o ud r u o ud )) + 4 a\ u u i dd / 7 (i«« ud r M i dd ) + 4 a™ d / 7 (o ud i" M i un o Md ) 

+ o;^ d0 " d (2 7 5 (0 ud l uu u<i l dd ) + 4/ 6 (l UM U(i ud l cM )) + a^ ^ (2 l 5 (l uu ud ud ud ) 

+ 4 j 6 (i««o" d O" d O" d ) + 4 / 7 (i« M o" d O Mc( Md ) + 4 j 7 (o Md l M "0 Mc( Mc! )) 

al dd ° ud = A + 2ari 4 (0 ud l dd l uu ) + 2 a\ M I 3 (l dd ud l dd ) + af "(2 / 3 (l dd u<i u<i ) + A l 4 (l dd ud ud ) 

+ 2/ 4 (0 MC W d )) + 4af "i dd / 7 (i^o^i^i««) + a i llu ° lld (2/ 5 (0 M(i 1^0 Md l nu ) 

+ 4/ 6 (l cW ud ud l uu )) + 4a 1 2 dd ° ud I 7 ( y ud l dd l dd ud ) + af d ° ud (2 l 5 (l dd ud ud ud ) 

+ 4 j 6 (i dd o Mc( Md O Mc( ) + 4 y 7 (i cW o u<i O U(i O u<i ) + 4 / 7 (o u<i l cM O u<i O u<i )) 

a°" d °" d = A + ^^(^(O^O^r") + 2 / 4 (0 Md Md l M ")) + ^(^(O^O^l^) + 2 / 4 (0 ud u<i l dd )) 

+ af i (2/ 3 (0"W) +4/4(0^0^0^)) + al uuldd (2I 5 (0 ud ud r u l M ) + 2 l 6 (0 ud ud l uu l dd ) 

+ 2 / 7 (o ud O ud l' i<i l uu ) + 2 j 7 (o Md O Md l M «l dd )) + al uu ° ud (2 / 6 (o Md O Md l M «0" d ) + 2 / 7 (o ud O" d l" M Md ) 

+ 2 j 7 (o« d o« rf o« rf i««)) + a 1 2 dd ° ud ( y 2 I 6 (G ud wi a ud l dd ) + 2 j 7 ( « W^) + 2 ^(o^O^l^O^)) 

+ af d01 " 1 (2 / 5 (o« rf O ud O ud O U(i ) + 2 / 6 (o« rf O ud O Md O Md ) + 4 / 7 (0 ud ud ud ud )) 

I ^jmijiiiiQiiii |2 j ^QudQud-^uuQud-^dd'j _|_ 2 / ^QudQuei^uit^ciciQMC^ -\- 2 I (Q u dQ u d-Qud-^uu 1^)) 

r 1 " 011 " = A + 2a\ uu I 9 (r u ud l dd r u ) + 2a\ dd I 9 (l dd ud r u l dd ) + a° ud (4 l 9 (l uu l dd ud ud ) 

+ 2 / 9 (i««o ud l dd O Md ) + 2 ^(l^O^lO™*)) + 4 a*"" 1 " J 10 (l W4 ,ld l <w l ,iU l <w ) 

+ 4 al uu ° ud l w (l dd l uu ud ud l uu ) + 4 af^'/oifl^O^O^l^) 

+ a f d ^ u \AI w {r u l dd G ad G ad G ad ) + 4/ 10 (i^i««o Md O Md O Md ) + 4/ 10 (l M "0 Md 1^0 Md M(i )) 
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f^S/j=_600 : 

a {" = A + 8 a\ ss h(l ss l ss ) + 12 a 1 ^ I 2 (l ss l ss l ss ) 

nf l '~ = A + af(4/ 3 (l s n ss l ss ) + 8/ 4 (i s n s n ss )) + al 33l33 (Ai 5 (i ss i ss i ss i ss ) 
+ 8 i 6 (i ss i ss i ss i ss ) + 16 / 7 (i s n ss i s n ss )) + 8 al 33133133 i 8 (i ss i ss i ss i ss i ss ) 

a i«i«i" = A . I2 a i ss / 9 (l ss l ss l ss l ss ) + 24 af 1 ' 3 1 10 (l ss l ss l ss l ss l ss ) 

AA SIJ= _ 200 : 

af" = A + AaT d h(l uu O ud ) + A af h(l uu O us ) + 2 af° ud l 2 (l uu ud ud ) + 8 af° us l 2 (l uu ud us ) 

+ 2aT 3 ° US h(l uu O us O us ) 

af = A + 4< d / 1 (l <i<i (r i ) + AaT 3 h(l dd O ds ) + 2 a^" I 2 (l dd O ud O ud ) + 8 af^^flW) 

+ 2af° d 7 2 (l <i<i O ds O ds ) 

a f = A + 4af7 1 (l ss O us ) + Aafh(l ss O ds ) + 2 af S °"7 2 (1 SS MS US ) + 8 af s ° d 7 2 (l ss O MS O ds ) 

+ 2 af ° d 7 2 (l ss 0^) 

af = A + ar/^O^r") + af^O^l**) + 2 af ^((T^) + 2 af 7 1 (O u<i O us ) 

+ 2 af /^(T^ 5 ) + ^^(O^tT*) + 2 af ^^((TW*) + 2 af d ° ds / 2 (0 u<i u<i ds ) 

+ 2 af s ° d 7 2 (0 nd MS ds ) 

af 3 = A + af7i(0 us l uu ) + ar/i(O us l ss ) + 2 af 7i (0 us ud ) + 2 af h(0 us us ) 

+ 2 afh(0 us ds ) + 2 af° us I 2 (0 us ud us ) + 2 af d ° ds J 2 (0 MS M( V s ) + af S °"7 2 (0 US US US ) 

+ 2 af S0d 7 2 (0 MS 0" s ds ) 

af = A + a{ dd h(O ds l dd ) + a{ 33 h(O ds l ss ) + 2 af h(0 ds ud ) + 2 af h( ds us ) 

+ 2 af 7 1 (0 ds ds ) + 2 af° U3 I 2 (0 ds ud us ) + 2 af° ds I 2 (0 ds ud ds ) + 2 af " 3 1 2 (0 ds us ds ) 
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+ af° ds I 2 (O ds O ds O ds ) 



= x + a i u "/ 3 ( ^o^r«) + a { l 3 {0 ud ud l dd ) + 2 a? I 3 {0 ud ud ud ) + 4 a?" 7 4 (0 U W s ) 

+ 4 af 7 4 (0" W s ) + 4 af d °"7 7 (0" d 0" d 0" d 0" s ) + 4 a°" "^^(O^O^O^O^) 

+ 2 s ° u 7 6 (0 ud 0" d us us ) + 4 af s ° ds (7 5 (0 ud ud us ds ) + ^(O^O^O^O*)) 

+ 2 ° d 7 6 (0 Md Md ds ds ) + 4 al ud ° us ° ds I 8 (0 ud ud us ud ds ) 

a ° 2 ud ° us = A + ar^(0 ud us r M ) + a° ud (/ 3 (0 ud us ud ) + 7 4 (0 us ud ud )) + aT° (h(0 ud us us ) 

+ h{0 ud us us )) + a°" 3 (I 3 {0 ud us ds ) + h{0 ud us ds ) + I 4 {0 us ud ds )) 

+ af d ° U3 \l 6 {0 ud us us ud ) + I 7 (0 ud us us ud ) + I 7 (0 us ud ud us )) + af dods (/ 5 (0 Md O MS U(i O ds ) 

+ I 6 (0 ud us ds ud ) + 7 7 (0 ud us M( V s ) + 7 7 (0 Mcf MS ds Md )) + s ° ds (7 5 (0 MS Md MS ds ) 

+ 7 6 (0 ud us us ds ) + 7 7 (0 us ud us ds ) + 7 7 (0" s 0" d ds us )) + af° ds I 6 (0 ud us ds ds ) 

af ¥ ' = A + a} dd /3(0 ud ds l dd ) + af d (/ 3 (0 ud ds ud ) + 7 4 (0 ds ud ud )) + a?" * (I 3 (0 ud ds us ) 

+ I 4 {0 ud ds us ) + 7 4 (0 ds ud us )) + af {I 3 {0 ud ds ds ) + 7 4 (0 ud ds ds )) 

+ a f d ° us (I 5 (0 ud ds ud us ) + I 6 {0 ud ds us ud ) + I 7 (0 ud ds ud us ) + I 7 (0 ud ds us ud )) 

+ af d ° d3 (I 6 (0 ud ds ds ud ) + 77(0^0^0^0^) + 7 7 (0 ds ud ud ds )) + a^ ^ h{0 ud ds us us ) 

+ al U3 ° dS (h(0 ds ud us ds ) + I 6 (0 ud ds us ds ) + 7 7 (0 ds ud us ds ) + 7 7 (0 ds ud ds us )) 

+ a 3 ,Ud ° US ° dS (I 8 (0 Ud dS ds ud US ) + I^dndsnusndsnud^ + j^udndsndsnusnud^ 

a „ s = x + Q ,i^/ 3 (0 US US 1 UU ) + a\" I 3 (0 us us r s ) + 4 d 7 4 (0" s MS ud ) + 2 a°" S 7 3 (0 US US 0" S ) 

+ 4 a?" 3 7 4 (0 us us ds ) + 2 af d ° ud I 7 (0 us us ud ud ) + 4 af dolls 7 7 (0 MS MS MS U(i ) 

+ 4a°" d0ds (7 5 (0" s us Md ds ) + I 6 {0 us us ud ds )) + 4 s ° ds 7 7 (0 MS MS MS ds ) 

+ 2 a° dS0ds 7 6 (0 MS MS ds ds ) + 2 a°" d ^*°* J 8 (0™0 M ,ld tta d8 ) 
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= A + ^"^(O^O*!") + af (J 3 (0 us ds ud ) + I 4 (0 us ds ud ) + h(0 ds us ud )) 

+ a°" 3 (I 3 (0 us ds us ) + h(0 ds us us )) + af (I 3 (0 us ds ds ) + h(0 us ds ds )) 

+ a° ud0Ud I 6 {0 us ds ud ud ) + a° ud ° us {I 5 {0 us ds ud us ) + I 6 (0 us ds ud us ) + I 7 {0 us ds ud us ) 

+ I 7 (0 us ds us ud )) + af d ° ds (I 5 (0 ds us ud ds ) + I 6 (0 ds us ud ds ) + I 7 (0 ds us ud ds ) 

+ I 7 (0 ds us ds ud )) + af s ° dS (h(0 us ds ds us ) + I 7 (0 us ds ds us ) + I 7 (0 ds us us ds )) 

+ afwv* (I 8 (0 us ds ds ud us ) + I 8 (0 us ds ud ds us ) + I 8 (0 us ds ds us ud )) 

= A + a{ dd I 3 (0 ds ds l dd ) + a{ 3a I 3 (0 ds ds r s ) + 4 a?" d h(0 ds ds ud ) + 4 a?" 7 4 (0 ds ds us ) 

+ 2<7 3 (0 ds ds ds ) + 2af d °" d / r (0 ,fa <fa u,i uc! ) + {h{0 ds ds ud us ) 

+ I 6 (0 ds ds ud us )) + 4 a° ud ° ds I 7 (0 ds ds ds ud ) + 2 a°" s °"7 6 (0 ds ds us us ) 

+ 4 s ° d 7 7 (0 ds ds 0^0 MS ) + 2 a 3 ,ud ° us ° ds I 8 {0 ds ds ud ds us ) 

= A + 7 9 (0 ud us ds r u ) + af d / 9 (0 ud ds us l dd ) + af7 9 (0 us ds ud l ss ) 

+ a? d (I 9 (0 ud us ds ud ) + I 9 (0 ud ds us ud ) + I 9 (0 us ds ud ud )) + a?" s (/ 9 (0 ud us ds us ) 

+ I 9 (0 ud ds us us ) + I 9 (0 us ds ud us )) + af s (I 9 (0 ud us ds ds ) + J 9 (0" d ds us ds ) 

+ I 9 (0 us ds ud ds )) + o;° udoud /io(0 us ud ds ud ud ) + a°" d ° us (I 10 (0 ud ds us ud us ) 

+ I w {0 us ud ds ud us ) + I w {0 us ds ud ud us ) + I 10 (0 ds us ud ud us )) 

+ af d ° ds (7 10 (0 ud us ds ud ds ) + 7 10 (0 ds ud us ud ds ) + 7 10 (0 ds us ud ud ds ) 

+ I w {0 us ds ud ud ds )) + a° us ° US ho{O ud O us O ds O us O us ) + af s ° ds (I 1Q (0 ud ds us us ds ) 

+ I w {0 ds us ud us ds ) + I w {0 us ud ds us ds ) + I w (0 ds ud us us ds )) 

+ af° ds I 10 (0 ud ds us ds ds ) 
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i '5/J=-3i2 : 

af = \ + 2afh(l uu ud )+6afl 1 (r u us ) 

af = \ + 2a 1 fh(l ss l ss ) + 4afh(l ss O us ) + 2 af h(l ss O ds ) 

af = A + a{ uu h(O ud r u ) + afh(O ud O ud ) + 3 af * h(O ud O us ) + 3 af h{O ud O ds ) 

af = A + af7 1 (O us l uu ) + 2al as h(0 us l ss ) + d h(O us O ud ) + 3af h(O us O us ) + af 7i(O us O dfl ) 

+ 2 a^" ls 7 2 (O us l uu l ss ) + 2 al 3a ° ud I 2 (O us O ud l ss ) 

af = A + 2a{ as h(0 ds l ss ) + 2 af h(O ds O ud ) + 2 af h(O ds O us ) + 2 af h(O ds O ds ) 

+ Aaf QUd I 2 (<d ds Q ud l ss ) 

: uiss = A + 2 af h{l ss l uu l ss ) + 2 af 7 4 (l uu l ss O ud ) + 4 af 7 3 (1 UU 1 SS (H 

+ 4a^ s °"7 6 (l uu l ss O ud l ss ) 

, ss °" d = A + 4 af h(O ud l ss l uu ) + 2 af h(l ss O ud r s ) + af I 4 (O ud l ss O ud ) + 2 af h(l ss O ud O us ) 

+ 2 af l 3 (r s ud ds ) + 2 af n 7 6 (l ss U(i l 5s r u ) + 2 a^ s ° u 7 6 (l ss O ud l ss O ud ) 

+ 2af 13aoda I 8 ( y l ss O ud l ss O ds l uu ) 

l-o* = A + 2ai s 7 9 (O ds l ss l uu l ss ) + 2af 7 9 (l uu O ds l ss O ud ) + af s (2 / 9 (r u ds l ss us ) 

+ 4/ 9 (l uu l ss ds us )) + 2af7 9 (O ds l ss l uu O ds ) +4ar°" d /io(l uu O a!s l ss O ud l ss ) 
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Table I. S-wave dibaryon masses. Parameters of model: cutoff A = ff.O, gluon coupling 
constants g = 0.653 and g x = 0.292. Quark masses m U:d = 410 MeV and m s = 557 MeV. 
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Table II. S-wave dibaryon masses. Parameters of model: cutoff A = 11.0, gluon coupling 
constants g = 0.647 and g x = 0.325. Quark masses m U:d = 410 MeV and m s = 557 MeV. 
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Table III. Vertex functions and Ghew-Mandelstam coefficients. 
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